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Refraction of shear zones in granular materials 
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We study strain localization in slow shear flow focusing on layered granular materials. A heretofore 
unknown effect is presented here. We show that shear zones are refracted at material interfaces in 
analogy with refraction of light beams in optics. This phenomenon can be obtained as a consequence 
of a recent variational model of shear zones. The predictions of the model are tested and confirmed 
by 3D discrete element simulations. We found that shear zones follow Snell's law of light refraction. 
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The motion of granular materials, such as sand, is dif- 
ficult to predict, as they behave neither like elastic solids 
nor like normal fluids. When they yield under stress and 
start flowing, the material consists of large, almost solid 
parts, and the relative motion is confined to narrow re- 
gions between them, called shear zones or shear bands 
flL S 9 • Shear zones represent material failure and there- 
fore play a crucial role in geophysics (geological faults), 
in engineering (building foundations) and in industrial 
processes. It is an important problem of these fields to 
predict where the failure takes place. This is a difficult 
task for those shear zones which arise in the bulk of the 
material far from the confining walls @, 0> H> B S 9 • 

Nature often stratifies granular media, i.e. different 
materials are deposited in distinct layers. We study here 
the question what effect such inhomogeneities have on 
shear zones. Based on a recent theory and computer 
simulations we show that shear zones are refracted when 
they pass through the interface between different granu- 
lar layers. Moreover, the angle of refraction obeys a law 
analogous to Snell's law of light refraction which reveals 
an unexpected analogy between granular media and ge- 
ometric optics. The effect of refraction presented here 
influences the position of shear zones in various kinds 
of materials including powders, sand, soil and rock lay- 
ers providing implications on the behavior of geological 

faults [1G3. 

Recently the variational principle of minimum dissipa- 
tion was successfully applied @, Hi to describe the non- 
trivial shape of the shear zone in a homogeneous granu- 
lar material which is sheared in a modified Couette cell 
0j & Q ■ Minimum rate of energy dissipation is a widely 
applied selection principle [Tl| to determine which steady 
states are realized in nature. When applying this prin- 
ciple to the cylindrical geometry of the modified Cou- 
ette cell , not only the non-trivial position of the shear 
zone could be calculated without any fitting parameter, 
but new, closed shear zones we re p redicted which were 
discovered also in experiments [l2| and simulations @|. 
In the following, a so far unknown but measurable effect 
will be derived from the same principle. It concerns the 
shape of shear zones in layered granular materials. 

This phenomenon can be understood best, if the side 



effects of gravity and curved shearing are absent. There- 
fore a long cylindrical container is considered, which is 
cut along its axis into two halves (Fig. QJi). (In com- 
puter simulations periodic boundary conditions in the 
axial direction are applied.) The container is completely 
filled with two granular materials having different fric- 
tional properties. They are separated by a planar inter- 
face parallel to the axis, but at an angle to the cut of 
the container. The system is sheared quasi-statically Q 
by moving the two halfs of the cylinder wall slowly par- 
allel to the axis in opposite directions. This creates a 
shear zone starting and ending where the container wall 
is cut. In order to find the shape of the zone in between 
we apply the principle of minimum dissipation. This, 
in zero width approximation Q, leads to the following 
variational problem (explained below): 
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(1) 



Here the local rate of energy dissipation is integrated over 
the whole surface of the shear zone. The shear zone is 
regarded as being infinitely thin, outside the shear zone 
no deformation and no dissipation takes place. The lo- 
cal dissipation rate per unit area vp/i eS is obtained by 
the sliding velocity v between the two sides of the shear 
zone times the shear stress p/i eS . The parameters p and 
// denote the overall pressure in the system and the 
coarse grained effective friction coefficient Q. The ef- 
fective friction has different values for the two materials 
while v and p are taken constant for our setup [1 91 ] . The 
surface that minimizes the total rate of dissipation pro- 
vides the shape of the shear zone. This means that the 
system yields along the surface where it has the least 
resistance against the external shear. 

The system is translation invariant along the axis, 
therefore the above surface integral is reduced to a line 
integral. We have to find the path (cross section of the 
shear zone) which connects two fixed points and along 
which the integral of p cS is minimal. This problem has 
exactly the form of Fermat's principle of geometric op- 
tics where the effective friction plays the role of the index 
of refraction. It is known that Fermat's principle of the 
shortest travel time leads to refraction of light beams 
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FIG. 1: (color online) a, Schematic view showing the shear cell 
and the refraction of the shear zone. Dense granular material 
is sheared between two half-cylinders. The shear direction is 
indicated by the arrows, b, A snapshot of the simulation. The 
upper light-gray (yellow) material and lower dark-gray (blue) 
material have different frictional properties (blue beads have 
stronger friction). If the velocity of any bead drops below 
half of the external shear velocity, it is overpainted by dark- 
gray (red). Despite the fluctuations of single bead- velocities, 
a dark (red) zone appears in the bulk that separates the two 
moving parts of the system, c, The spatial distribution of 
the shear rate which shows the structure of the shear zone. 
Lighter colors represent stronger deformations (orange means 
the largest shear rate, followed by red and blue, while black 
means almost no deformation). The straight line indicates 
the material interface. 



therefore one can expect a similar phenomenon of refrac- 
tion also for shear zones. 

In order to check this prediction we carried out com- 
puter simulations of 10 frictional and hard spherical 
beads (Fig. [TJd). Based on a standard discrete element 
method (DEM) [HB3 

the motion of each bead is deter- 
mined by Newton's equation. This grain-level dynamics 
leads to a collective steady state flow where we measure 
the coarse-grained velocity field. The coarse-grained ve- 
locities are parallel to the cylinder axis. The local strain 
rate, which is calculated from the velocity field, corre- 
sponds to pure shear deformation (without volumetric 
strain). The spatial distribution of the magnitude of the 
shear rate is shown in Fig. [TJ;. As predicted, the arising 
shear zone departs from the cut plane and its direction is 
changed significantly at the material interface (Fig. [TJd, 
c). Note, that the shear zone we obtained is quite wide 
relative to the system size. The relative width becomes 
smaller for larger systems [3| which, however, would de- 
mand also a computational effort out of reach. Still, the 
effect of refraction is clearly shown. 

The analogy with light refraction provides quantita- 
tive predictions that can be measured in numerical and 
real experiments. Fermat's principle indicates that shear 
zones have to follow Snell's law of light refraction which 
will be tested here by DEM simulations. In order to do 
so we characterize the behavior of the shear zone by two 
angles (f>\ and 4>2 which are analogous to the incident and 
refractive angles used in geometric optics. These are de- 
fined with the help of the shear deformation (Fig. [TJ;) 
as follows: first, the refraction point is identified as the 
point along the interface where the local shear rate has 
its maximum. Then the refraction point is connected 
to the starting and ending points of the shear zone by 
straight lines. Thus the direction of the shear zone on 
each side of the interface is determined. These directions 
provide the angles <p\ and <p2 with respect to the normal 
of the interface. Another angle, which is an important 
input parameter of the shear test, is the tilt angle of the 
interface a between the normal of the cut plane and the 
interface (Fig. [JJ;). For various values of a we measured 
</>i and 4>2- If the interface is perpendicular to the cut 
plane {a = 0) no refraction was observed, however, as 
a is increased the refraction of the shear zone becomes 
more and more pronounced. 

If we apply Snells law to our case it indicates that, 
no matter how the angles </>i and 02 are changed, the 
ratio of sin 4>\ and sin 02 has to be independent of the 
tilt angle a. Furthermore it is expected that this ratio 
can be expressed by the effective frictions: 

sin 0i _ A^*2 (0 -s 

• X pff v > 

sin <p2 Mi 

First we report the influence of a observed in simula- 
tions. We used materials that contained hard spheres of 
radii distributed uniformly between i? m i n and 1.3i? m j n . 
This polydispersity was needed to avoid shear induced 
crystallization 15]. In order to control the frictional 
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FIG. 2: (color online) 4>i (a), 02 (b) and the relative index of 
refraction (c) measured for shear zones in computer simula- 
tions. The horizontal axis show the tilt angle of the interface 
between the two materials. Data points denoted by triangle, 
square and circle stand for systems A, B and C, respectively. 
The pairs of microscopic friction that are used in these sys- 
tems are indicated in the figure. The straight lines show the 
theoretically predicted values (the upper line for system A 
and B, the lower one for system C). 



properties of the materials we varied the value of the mi- 
croscopic friction /i which is the friction coefficient at the 
particle-particle contacts, ft is important to note that 
the microscopic friction \i and the effective friction /x eff 
are not the same although they are closely related as we 
will see later. We simulated three systems called A, B 
and C. The total number of particles was 1 00000 for sys- 
tem A and 50000 for systems B and C. For each case the 
shear cell was filled with two materials with microscopic 
frictions [i\ and fi 2 - For systems A and B [i\ = and 
Hi = 0.5 while for system C fii = 0.1 and fi 2 = 0.5. For 
all systems the radius of the container was approximately 
65 R 

mm ■ 

The influence of the tilt angle of the interface is pre- 
sented in Fig. The computer simulations show that <f>\ 
and 4>2 vary strongly but sin cf)\ / sin <p2 remains constant 
for a wide range of the angle a. This ratio seems to de- 
pend only on the materials in which the shear zone was 
created, in full agreement with the theoretical consider- 
ations. 

Next we deal with the role of the effective friction. It 
is known that if a dense noncohesive granular material is 
sheared then after a short transient [161 ] it reaches a well 
defined resistance against shear deformation provided the 
flow is slow enough (quasistatic shear). Thus a mate- 




FIG. 3: Effective friction versus microscopic friction. The 
open circles correspond to the materials used in systems A, 
B and C. 



rial parameter, the effective friction coefficient, can be 
defined which turns out to be independent of the prepa- 
ration history, the actual shear rate and the confining 
pressure [!, [Tj]. The effective friction [i cS is given by 
the shear stress divided by the normal stress, both mea- 
sured in the plane along which the shear deformation 
takes place. 

The quantity /i eS is not a microscopic input parame- 
ter of the computer simulation, thus the question arises 
what values fj, has for the materials that appeared in 
systems A, B and C. This can be deduced from simu- 
lations which are independent of the previous refraction 
tests: we put the same materials into a 3D rectangular 
box under plane shear 0] measure the components of the 
stress tensor [l8j and calculate the value of /i Gff . fn this 
way we achieve a calibration curve that provides a one 
to one correspondence between microscopic and effective 
friction, see Fig. [3] 

With that we arrived to a point where we can make 
quantitative statements about the extent of the refraction 
based on Snell's law. E.g. for system A, with help of 
the prescribed values of \i\ and /X2, the ratio ^2/^1 
can be calculated and compared to the data sin <f>\ / sin <j}% 
recorded in the refraction test. This is done in Fig. [2] for 
systems A, B and C. A surprisingly good agreement is 
found between the predicted and measured values which 
holds for all systems and for various tilt angles of the 
interface. This is our main result. 

Dealing with shear localization in layered granular ma- 
terials we investigated the predictions of a recent varia- 
tional model and compared them to computer simula- 
tions. The results of the present work convey two mes- 
sages. First, we gave verification of the variational ap- 
proach to granular shear flows. We tested the model in 
a new situation and, according to the numerical data, it 
gave an excellent description of the behavior of the shear 
zones. Second, an interesting analogy between granular 
flow and geometric optics is revealed. We showed that 
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shear zones are refracted at material interfaces similarly 
to light beams. The phenomenon presented here should 
be accessible by experiments. 
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